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ABSTRACT 

We present spectroscopy of binary quasar candidates selected from Data Release 4 of the Sloan 
Digital Sky Survey (SDSS DR4) using Kernel Density Estimation (KDE). We present 27 new sets of 
observations, 10 of which are binary quasars, roughly doubling the number of known g < 21 binaries 
with component separations of 3" < A8 < 6". Only 3 of 49 spectroscopically identified objects are 
non-quasars, confirming that the quasar selection efficiency of the KDE technique is ~ 95%. Several 
of our observed binaries are wide-separation lens candidates that merit additional higher-resolution 
observations. One interesting pair may be an M star binary, or an M star-binary quasar superposition. 
Our candidates are initially selected by UV-excess (u — g < 1), but are otherwise selected irrespective 
of the relative colors of the quasar pair, and we thus use them to suggest optimal color similarity 
and photometric redshift approaches for targeting binary quasars, or projected quasar pairs. From a 
sample that is complete on proper scales of 23.7 < R pr0 p < 29.7 h^ 1 kpc, we determine the projected 
quasar correlation function to be W p = 24.0±iq;§, which is 2a lower than recent estimates. We argue 
that our low W p estimates may indicate redshift evolution in the quasar correlation function from 
z ~ 1.9 to z ~ 1.4 on scales of R pr0 p ~ 25 h^ 1 kpc. The size of this evolution broadly tracks 
quasar clustering on larger scales, consistent with merger-driven models of quasar origin. Although 
our sample alone is insufficient to detect evolution in quasar clustering on small scales, an i-selected 
DR6 KDE quasar catalog, which will contain several hundred z < 5 binary quasars, could easily 
constrain any clustering evolution at R prop ~ 25 h^ 1 kpc. 

Subject headings: cosmology: observations — large-scale structure of universe — quasars: general — 
surveys 



1. INTRODUCTION 

Cosmologically, quasars can now be explained as one 
spectacular stage of an evolutionary process, initiated by 
gas-rich galaxy mergers, th at ultimately helps redden e l- 
liptical galaxies (see, e.g., iHopkins et al.l [20061 [2007a). 
That quasar activity might trace the early stages of 
merger-driven galaxy evolution makes quasar observa- 
tions an essential ingredient in constraining galaxy for- 
mation scenarios. On the other hand, less luminous Ac- 
tive Galactic Nuclei (AGN), particularly at low redshift 
(z < 1), may be bett er explained by less violen t fuel- 
ing mechanisms (e.g.. [Hopkins fc Hernquistll2006l ) than 
the mergers that drive quasars' optical intensity. Given 
that orientation can complementarily explain many dif- 
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ferences be tween the AGN zoo (e.g., lAntonuccl 119931 : 
new AGN constraints over a broad range of 
luminosity, particularly across z ~ 1, may be key to de- 
termining which elements of quasar behavior are mainly 
structural, and which are mainly evolutionary. 

If quasars are associated with galaxy mergers, obser- 
vations of binary quasars with proper separations com- 
parable to the scale of small galaxy groups should of- 
fer interesting constraints. It has generally become ac- 
cepted that most quasar pairs at similar redshifts that 
have image s eparations > 3", are binary quasars rather 
than lenses dPhinnev fc Blandfordl Il986t IBahcall et al.l 
[19861: iKochanek et all 119991 : iRusinl l200a iQguril 120061) . 
AGN activity can be exacerbated by tidal forces in 
galax y mergers ([Barnes fc Hernquistlll996l : IBahcall et al.l 
119971) and it has long been argued that this might 
explain an excess of binary quasars flbiorgo vski 1991; 
IKochanek et al . 1999; Mort lock etalll 19991) and physica l 
triples (|Diorgovski et al.ll2007[ ). IHopkins et ail (|2007bf ) 
suggest instead that a binary quasar excess simply re- 
flects the increased probability for mergers to occur in 
regions that are overdense on small scales. If quasars 
form in merger s they will thus be nat urally more biased 
at small scales. IHopkins et al.l ([2007b) further argue that 
orbits for which quasar activity might be exacerbated in 
both of two merging galaxies are prohibitively rare, even 
if few s uch events are need ed to explain a binary quasar 
excess ([Myers et al.ll2007bt henceforth M07b). 

The Sloan D igital Sky Survey (henceforth SDSS; 
lYork et al. 2000) has renew ed interest in binary quasars, 
and iHennawi et alj ([20061 : henceforth H06) have used 
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SPSS data to categorically conf i rm ea rlier evidence (e.g., 
iDiorgovskil 119911 : iHewett et all 11998ft that quasar clus- 
tering is enhanced on comoving scales < 100 h^ 1 kpc 
(proper scales of < 40 ft. -1 kpc at z ~ 1.5). Binary 
quasars are scarce, with perhaps fewer than five hundred 
in the entire sky to g < 21, redshifts of z < 2.5, and 
comoving separations < 100 hr 1 kpc (given the sample 
size in M07b). Deep, wide imaging, such as from the 
SDSS, is thus key to testing predictions of the nature of 
binary quasars; for instance, by studying quasar cluster- 
ing as a function of redshift or luminosity. If enhanced 
small-scale quasar clustering is due to the enhanced bias 
of major galaxy mergers, rather than tidal forces exac- 
erbating quasar activity, then there should be no red- 
shift evolution in the relative bias of quasar clustering 
on large and small scales. Further, some merger-driven 
models predict stronger small-scale clustering for quasars 
than for lower-lu minosity AGNs. due to different fueling 
mechanisms (e.g., Ho pkins et al.ll200"7bt ). 

Binary quasar samples are typically assembled from 
pairs that a re initially targeted a s possible gravitational 
lenses (e.g.. iMortlock et"aLl il999: H06). Because of this, 
the two members of most known binary quasars have 
similar colors. Although color similarity may be optimal 
in detecting l enses, scatter in the quasar color-redshift 
relation (e.g., Richard s et al.l l2001h dictates that strict 
color similarity cannot select all binary quasars. As bi- 
nary quasars are useful in testing merger-driven models 
of quasar activity, it is disconcerting that color similar- 
ity cuts might discard particularly informative binaries, 
such as any that are being exacerbated by tidal forces in 
merging galaxies. Further, as binary quasars are scarce, 
relaxing color criterion might help provide enough bina- 
ries to study their redshift evolution. 

We have now performed spectroscopy of a sam- 
ple drawn from a complete, UVX (u — g < 
1) set of SDSS Data Releas e 4 (henceforth DR4; 
lAdelman-McCarthv etldl 120060 binary quasar candi- 
dates (see Table 1 of M07b). The sample, ~ 45% of which 
has now been identified, is photometrically selected us- 
ing the Kernel Densi ty Estimation (KDE) technique of 
IRichards et all (|200l . Our aim is to compile an exten- 
sive, homogeneous set of binary quasars, with proper sep- 
arations < 40 h~ x kpc, mainly to study quasar cluster- 
ing on small scales, ultimately as a function of redshift. 
Most of our binary candidates are unlikel y to be observed 
in th e main SDSS quasar survey (e.g., ISchneider et al.l 
120071) . because our observations probe fainter, and be- 
cause SDSS fibers (in single tiles) cannot be placed closer 
than 55". Our approach differs from previous stud- 
ies of binary quasars. In particular, our main goal is 
to study binary quasars, not gr avitational lenses (un- 
like the samples com piled in, e.g.. iKochanek et"aLlll999l ; 
IMortlock etaLlll999t H06). Thus, excepting our initial 
UVX cut, our binaries are the first sample selected re- 
gardless of the relative colors of the component quasars. 
This allows an investigation of whether color similarity 
can be used to optimally target binary quasars ($3|). 

In we discuss our initial results in compiling a ho- 
mogeneous, spectroscopic binary quasar sample from the 
SDSS, and report 27 sets of new observations of binary 
quasar candidates. In f]4]we use a subset of these observa- 
tions to present the first analysis of quasar clustering on 



small scales using a complete spectroscopic sample of bi- 
nary quasars. All our binaries are selected from the DR4 
KDE candidates of M07b, and are therefore a straightfor- 
ward subset of SDSS DR4, making our selection function 
very simple. We correct al l magnitudes for Galactic ex- 
tinctio n using the maps of ISchlegel. Finkbeiner fc Davisl 
(1998), unless otherwise noted. We adop t fl m = 0.27, 
fl A = 0.73, consistent with WMAP3 (|Spergel et all 
l200l . and h = iJ /100km s" 1 Mpc" 1 = 0.7. We de- 
note transverse proper (comoving) scales as R prop (R) 
and radial proper (comoving) scales as s prop (s). 

2. DATA 

2.1. Observations 

2.1.1. Candidate Selection 

Our candidate binary quasars are photometric objects 
in SDSS D R4 that are classified as quasars by the KDE 
technique (Ri chards et al.ll2004T) . have u — g < 1, g < 21, 
and are within 3" to 6" of ano ther such object . The SDSS 
ugriz filters are described in iFukugita et all (|1996l ). As 
in M07b, we inspect our candidates and discard any 
that are clearly not quasars (close KDE candidates are 
occasionally misclassified H II regions in low-rcdshift 
galaxies). We restrict our spectroscopy, and analysis, to 
quasar pairs with angular separations of Ad > 3", so that 
components of a pair appear clearly separated in SDSS 
imaging (e.g., H06). Although scales that are comparable 
to small galaxy groups or smaller are most useful in test- 
ing merger models of quasar activity our upper limit of 
A(9 < 6" is somewhat arbitrary, providing a reasonable 
number of candidates for a small spectroscopic study. 
Table 1 of M07b lists our 98 DR4 candidate binaries (and 
a further 13 candidates with A9 < 3"). At g < 21 the 
KDE technique, coupled with a UVX (u — g < 1) cut, 
selects quasars with 95% efficiency, and is over 95% com- 
plete for redshifts of 0.2 ^ z < 2.4 (e.g.. IRichards et all 
l200l iMvers et al.l 120061 |2007a). Our survey goal is an 
efficiently classified, statistical, sample of binary quasars, 
which are UVX but are otherwise selected irrespective of 
the relative colors of the components of the binary. At 
z < 2.5, our selection should thus only bias our sam- 
ple against those binaries where one or both component 
quasars is reddened beyond u — g = 1 . 

2.1.2. Spectroscopy 

Spectroscopy of our DR4 KDE binary quasar candi- 
dates was obtained with the R-C Spectrograph on the 
Mayall 4-m, over 5 nights (UT 2007 February 22-26) 
at Kitt Peak National Observatory. We used a 1.5" by 
98" long-slit set at the position angle of the candidate bi- 
nary, allowing both components to be simultaneously ob- 
served. The KPC-10A grating and WG360 blocking filter 
yielded a resolution of ~ 5A and a wavelength coverage 
of ~ 3800-7500A. The most useful observations were ob- 
tained on February 25 and 26, due to cloud and wind 
on other nights. Over February 25 and 26, the seeing 
was < 1.5", allowing even our closest binary candidates 
(^3") to be spatially separated. The survey goals were 
a positive identification and redshift for each DR4 KDE 
binary candidate. This typically required a 15 minute 
exposure when the faintest member of the candidate bi- 
nary was at g ~ 19.5 and three 20 minute exposures 
(or longer) for g ~ 21, although the > 70% illuminated 
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Moon on February 25 and 26 typically prevented our 
g ~ 21 candidates from being spectroscopically identi- 
fied. 

Spectra were reduced at the telescope, using IRAF 9 . 
Exposures ceased once a binary candidate could be iden- 
tified as; (1) containing one star or galaxy; or (2) con- 
taining two quasars with established redshifts. We esti- 
mate redshifts using the rest frame emiss i on wa velengths 
listed in Table 2 of iVanden Berk et all (|2001h . Based 
on deriving redshifts, where possible, from several differ- 
ent lines in each quasar's spectrum, our typical error is 
Az ~ 0.0031. A similar approach suggests a mean ve- 
locity precision of ~ 370 km s^ 1 , with little redshift de- 
pendence. We note that this is likely an overestimate, as 
velocity differences between quasars are usually found to 
be more precise when cross-correlating t he full spectra, 
rather than measuring line shift s (e.g., iTonrv &; Davisl 
1979; Diorgo vskTfc SpinradlH98l . 

Tables [TaHldl detail our new observations. Table [Ta] 
lists objects for which we obtained an identification for 
only one member of the candidate binary. Table llbl lists 
confirmed binary quasars. Following H06, we classify 
quasar pairs with a line-of-sight velocity difference of 
|Au||| < 2000 km s" 1 (s prop < 9.5 hr 1 Mpc at z = 1.4) 
as a binary. Table [Tc] lists DR4 KDE binary candidates 
that are actually "projected" quasar pairs (with compo- 
nents that lie at different redshifts), or pairs in which 
one object is not a quasar. When both redshifts listed 
in Table [Tc] have some uncertainty, a binary quasar in- 
terpretation can still be ruled out, based on strong lines 
observed in both quasars but at discrepant wavelengths. 
Pairs that are harder to definitively identify are listed in 
Table ED (see g2JL2])- 

In all, we spectroscopically identified 49 objects, 
mainly in areas with Galactic absorption A g < 0.17. Of 
these identified objects, 44 are both members of 22 candi- 
date quasar pairs, and 5 are objects from pairs for which 
we identified only one member. Of the 49 identified ob- 
jects 46 are quasars, confirmin g the KDE techniq ue is 
~ 95% efficient for A g < 0.21 (jMvers et al.ll2007an . Of 
the 22 candidate quasar pairs for which we identified both 
components, 3 are quasar-non-quasar pairs, 9 are pro- 
jected quasar pairs (i.e. at disjoint redshifts), and 10 are 
binary quasars. Several of the 10 binary quasars could, 
in fact, be previously unrecorded lenses; see N2.2.1j) . In 
Table [2] we list previously known quasar pairs that also 
meet the criteria to be included in our DR4 KDE binary 
candidate sample. Previously known candidates include 
12 non-binaries, 9 binaries, and 2 lenses. Approximately 
half of the 98 (3" < A9 < 6") DR4 KDE candidates have 
now been spectroscopically identified (see Tabled]). With 
the caveat that bright objects may have been observed 
first, ~ 42% of the candidate pairs are binary quasars, 
and only ~ 16% of the pairs contain a non-quasar. 

2.2. Interesting Spectroscopic Pairs 
2.2.1. Potential Lenses 

Five pairs in Table E3 have sufficiently similar spec- 
tra, at our ~ 5A resolution, that they might be a lensed 

9 Distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 



quasar rather than a binary. As lenses with image sep- 
arations in the r ange 3" < A9 < 6" are rare (e.g., 
Ilnada et~aT1l2007t ). particularly f or z < 2, we interpre t 
these objects as binaries (see also iKochanek et al"1ll999|h 
although they certainly merit higher-resolution spec- 
troscopy. A lensing interpretation is especially unlikely 
for the three possible lenses used in our clustering anal- 
ysis (j2J, which have ~ 5" separations and, in two cases, 
dissimilar colors. In Figure [JJ we display the spectra of 
our most likely lens candidates. SDSSJ1158+1235A and 
B, in particular, have almost identical spectra at our res- 
olution. 

2.2.2. Notes on Ambiguous Binaries 

Table [Tdl lists four candidates that we could not defini- 
tively identify. We conclude that two of these objects 
are binaries, for the following reasons, quoting all wave- 
lengths in the observed frame. 

SDSSJ093424.32+421130.8 and 
SDSSJ093424.11+421135.0 consist of a quasar at 
z = 1.339 and a featureless spectrum (after 4200s of 
exposure). SDSSJ093424.11+421135.0 is faint (observed 
g = 21.01), probably a star, and has no obvious emission 
near 4460A or 6550A, the principal emission lines used 
to identify SDSSJ093424.32+421 130.8. We therefore 
conservatively conclude that this is not a binary quasar. 

SDSSJ120727.09+140817.1 and 
SDSSJ120727.25+140820.3 both have ambiguous 
redshifts. The fainter object (observed g = 20.39) has 
broad emission at 4320A, and near 7840A at the red 
edge of our coverage. Although the brighter object 
(g = 20.27), has possible, low signal-to- noise ratio, 
emission near 4340A we discount it based on the 
4320A emission in the fainter object being strong, and 
that the objects are, obviously, observed under similar 
conditions. We tentatively conclude that this pair is not 
a binary quasar. 

SDSSJ1235+6836A,B (see Figured]) are an interesting 
pair with highly dissimilar colors. SDSSJ1235+6836B 
is apparently a quasar with significant broad emission 
near 3895A and 7050A and weaker emission near 4790A. 
SDSSJ1235+6836B is most likely a quasar, with probable 
broad emission near 3910A and possible emission near 
4820A, lying behind a classic M star observed at the red 
end of the spectrum. As < 3900A, is near the blue edge 
of our coverage, an alternative possibility is that one or 
both objects are just M stars with very strong, blended 
Call H and K emission. Assuming the proper motion of 
the M star(s) is moderate, this object may look more or 
less like an M star-quasar pair over time. 

SDSSJ1507+2903A,B have strong emission near 
5250A and 5215A, respectively. Neither spectrum has 
additional features over 3800-7500A, and so we assume 
that the emission is Mgll, placing both quasars at z ~ 
0.87. The ambiguity for this pair is that their redshifts 
imply \Av\\\ = 2100 km s" 1 . As a shift of Sz < 0.0005, 
far smaller than our typical precision, can bring these 
quasars within \Avn\ < 2000 km s" 1 , we identify this 
pair as a binary. 

3. COLOR SELECTION OF BINARY QUASARS 

To study relative color selection of binary quasars, we 
use the \ 2 color similarity statistic introduced by H06 
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where the subscripts 1 and 2 represent the components 
of a pair. The superscript i refers to flux (/) in the 5 
SDSS bands (ugriz). For asinh magnitudes (m) 
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where P = 2.5/ In 10 (|Pogsonll856f ). F = 3 630.78Jy, and 
h [u, g ,r y] = [1,4, 0.9, 1 . 2, 1.8 , 7.4] x 10" 10 ([Lupton et all 
1999: iStoughton et al.ll2002h . A quasar pair with more 
similar colors has a lower XT&ior- Iterative equations for 
calculating A in Equation [l] can be ill-conditioned for 
Xcoior ~ ^0, so, throughout this work, we numerically 
determine A by bisection. 

Binary quasars are often pairs rejected from grav- 
itational lens searches, and as such, the components 
of known binaries typically have very similar col- 
ors, a long-known example being SDSSJ1637+2636A,B 
(yL_ = 2.8; see Table E3 iSramek fc Weedmanl Il978l: 
iDjorgovski fc Spinradlll984T K Schemes designed to opti- 
mize binary quasar searches by selecting pairs with simi- 
lar colors, will, therefore, naturally reselect known binary 
quasars. Our DR4 KDE objects are simply all candidates 
with a high probability of being quasars and thus, after 
the initial homogeneous UVX cut, are selected irrespec- 
tive of the relative colors of the components of the pair. 
The UVX cut itself, at z < 2.5, should only bias our 
sample against those binaries with a component that is 
intrinsically dust-reddened beyond u — g — 1 . Our sam- 
ple should thus be useful in determining color similarity 
cuts to optimize binary quasar selection. However, some 
of the quasar pairs in Table [5] were selected by H06 to 
have Xcoior < 20; as we avoided reobserving these pairs, 
our data in Tables fTblTTdl may be biased to x\ ior > 20- 

In the upper-left panel of Figure [31 we demonstrate 
that the Xcoior values of the DR4 KDE binary candi- 
dates identified to date fairly represent the full sam- 
ple. We compare the cumulative fraction of the 45 
identified candidates (i.e., Table [3]) to the remaining 53 
(3" < AO < 6") candidates, as a function of Xcoior- A 
two-sample Kolmogorov-Smirnov test cannot distinguish 
the distributions, suggesting that the colors of the ob- 
served candidates fairly represent all candidates. The 
upper-right panel of Figure [3] compares the Xcoior cu ~ 
mulative probability for the 21 confirmed binary quasars 
(or lenses) and the 24 confirmed non-binaries (projected 
quasar pairs, star-quasar pairs, NELG-quasar pairs). 
The K-S test probability that these two distributions are 
drawn from the same underlying Xcoior distribution is 
~ 10%, suggesting that XjLi™, can indeed discriminate 
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binary quasars from non-binaries. 

As our candidates are selected without a Xcoior cu t> we 
can ask what Xcoior limit optimizes completeness (num- 
ber of binaries/21 total binaries) and efficiency (num- 
ber of binaries/45 observed candidates). To sample 
> 50% (> 66%) of binary quasars requires Xcoior ~ 10 
ixloior £ 20), efficient at 70% (55%), while xLor £ 70- 
100, contains 95% of all binaries (only missing the "am- 



biguous" SDSSJ1235+6836A,B) and remains 50% effi- 
cient. A xloior < 70 cut rejects all but one quasar-star 
pair, while retaining all quasar-quasar projections. 

As there is reasona ble scatter in the q uasar color- 
rcdshift relation (e.g.. iRichards et"aT1l2001[ ) full photo- 
metric redshift (henceforth photoz) information should 
better select binary quasars. To te st this, we consider the 
primary "CZR" photoz solution (|Weinstein et al.l 12004 
"Zphot range" in Tables 1-2) for each quasar in our sam- 
ple. We determine the overlap fraction of the primary 
photoz solutions of the two quasars in a pair, multiplying 
by the probability that the quasar occupies that primary 
peak. 

The completeness and efficiency of a binary quasar 
sample obtained by considering photoz overlap are plot- 
ted in the lower-left panel of Figure [3] Confirmed non- 
binaries typically have no photoz overlap. A proba- 
bility cut at > 3% overlap will return 90% of bina- 
ries and is 73% efficient. The two binaries that are 
missed are the "ambiguous" SDSSJ1235+6836A,B, and 
SDSSJ1637+2636A,B, the "A" component of which has 
a poorly behaved photoz solution. If one additionally ob- 
served all candidates that contained a quasar with a poor 
photoz (characterized by a probability of < 0.5), 95% of 
binary quasars would be observed at 69% efficiency. Of 
course, although cutting on photoz overlap is a more ef- 
ficient mechanism for selecting binaries, it does so at the 
expense of projected quasar pairs. A cut at < 3% over- 
lap could therefore be used to discard binary quasars in 
favor of projected pairs. 

In conclusion, if a survey's goal is to select both binary 
quasars and projected quasar pairs, a cut of Xcoior ~ 70 
(after first applying an efficient photometric quasar clas- 
sification technique such as KDE) will return 97% of bi- 
naries, lenses and projected pairs (multiplied by, e.g., 
the 95% completeness of the KDE technique itself), at 
93% efficiency. This is hardly surprising, and simply fur- 
ther confirms that the KDE technique i s ~ 95% efficient 
(jRichards et al.ll2004t Mvers et alj|2006ll2007al ). Cuts of 
Xcoior ~ 20, or stricter, may be necessary without prior 
KDE photometric classification but will miss > 33% of 
all binary quasars. Interestingly, 2 (~ 10%) of the DR4 
KDE binary quasars (and one lens!) have more dissimilar 
colors than any of the quasar-quasar projections, hinting 
that physical interactions may affect the colors of a few 
binary quasars. If a survey goal is observing only binary 
quasars, a cut of > 0.03 in the overlap of the two quasars' 
photozs is superior to Xcoior selection. We stress that this 
analysis applies only to quasars pre-selected using an ef- 
ficient photometric classification technique. 

4. PROJECTED QUASAR CLUSTERING AT 25 h^ 1 kpc 

With our new observations ( Table [TBI and fid)) . we have 
now identified all DR4 KDE binary quasar candidates 
with An < 0.17, g < 20.85 and 3.9" < A0 < 5.2". In 
Table 0] we compile the binaries that meet these crite- 
ria. The three possible lenses in this subsample are wide- 
separation (A9 > 4.5"), making them likely binaries. As 
our sample is complete over 3.9" < Ad < 5.2", it is also 
effectively complete for quasars in SDSS DR4 with sepa- 
rations on proper scales of 23.7 < Rp r0 p < 29.7 Yr 1 kpc 
over 1.0 < z < 2.1. 

We st udy quasar clusterin g using the DD / DR estima- 
tor (e.g., Sha nks et al .11 1983? ) for quasar-quasar (QQ) pair 
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counts compared to expected quasar-random (QR) pair 
counts 



jit QQ 



w„ = 



(QR) 



(3) 



Higher-order corrections fe.g. JLandv k. Sz alav 1993|), re- 
duce to Equation [3] for the small scales and large vol- 
umes we consider . We use small-number Poisson er- 
rors from iGehrelsl (|1986f h Poisson errors are valid on 
small scales where pa i r counts are inde pendent (e.g., 
iCroom fc Shajikill996t iMvers et al.ll2006h . 

As our clustering sample is a subset of all SDSS DR4 
photometric objects, our selection function is very sim- 
ple. We calculate (QR) in Equation [3] by constructing 
a catalog of random points with the same angular cov- 
erage as S PSS DR4, cor r ecting the SP SS data for mask 
holes, as in IMvers et "all l|2006tl2007al) . We further limit 
our random catalog to areas of the sky with Galactic ab- 
sorption A g < 0.17. We assign random points a redshift 
according to a fit to the normalized redshift distribution 
of (A„ < 0.17, q < 2 0.85) quasars in the PR1 catalog 
([Schneider et al.ll2003h . from which t he PR4 KPE quasar 
classification is trained. Figure 7 of IMvers et all ([2 006) 
is sim ilar to this redshift distribution, and IMvers et al.l 
(2006) argue that including additional quasars that over- 
lap the KPE color space minimally impacts this N(z) 
distribution. To represent our fit, we use a modified 
Gaussian 



dN = ftexp- 



-dz 



no- 
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and find a best fit with z = 1-4, o = 0.6, ft = 0.65 and 
n — 3 (see also IMvers et al.l [2007af) . We have repeated 
our analyses instead using a spline fit, and our results 
differ by < 2%. 

To determine (QR) in a given bin we use a random cat- 
alog with 1000 times as many points as the (A g < 0.17, 
g < 20.85) KPE PR4 photometric quasar catalog. We 
total all QR counts in the angular bin of interest, and 
normalize the result (divide by 1000). We then create 
a random catalog with points distributed according to 
Equation |4j and determine the fraction of pairs that 
would lie within |AuJ < 2000 km s _1 by Monte Carlo 
sampling to 0.1% precision. Multiplying this fraction by 
the normalized QR counts yields (QR). In a bin of 3.9" < 
A0 < 5.2" we expect 30.3 QR counts. Over our full N(z), 
we expect 0.0166 pairs with |Auii | < 2000 km s -1 . Thus 
(QR) = 0.503, compared to QQ — 16 for the 8 (non- 
unique) pairs in Table [H The implied projected corre- 
lation function averaged over 3.9" < A9 < 5.2" is thus 
W p = 30.8±Jf'Q. We note that, for all candidate pairs, 
QQ(3.9" < 8 < 5.2") = 40, yielding an angular corre- 
lation of w(6) = (40/30.3) - 1 = 0.320±g |g|, consistent 
with M07b. 

One of the pairs (SPSSJ1507+2903) in Table H has 
|Aii||| < 2100 km s . We reasonably include this pair 
in our analysis, given the precision of our redshift esti- 
mates. Instead rejecting SPSSJ1507+2903 and assum- 
ing QQ = 14 would lower our estimate of W p by ~ 10%. 
The |A«||| < 2000 km s _1 limit from H06 is intended 
to bracket possible shifts in quasar lines, and thus in- 
corporate possible errors on |Ai>|| | but is otherwise arbi- 



trary. Although our Monte Carlo sampling of |Au|| | < 
2000 km s _1 pairs does not model any error in |Avii|, 
relaxing this velocity window to |Av||| < 2100 km s _1 
would imply (QR) = 0.533 for QQ = 16, lowering our 
W p estimate by ~ 5%. As the change in W p implied 
by relaxing our velocity criterion is far smaller than our 
errors on W p , we proceed including SPSSJ1507+2903 in 
our analyses and maintaining \Avn\ < 2000 km s _1 in 
our Monte Carlo sampling. 

We apply our approach to different redshift ranges 
by weighting the number density of points in the ini- 
tial calculation of QR by the relative fraction of PR4 
KPE quasars obtained by integrating under Equation 01 
We can thus determine (QR) over 1.0 < z < 2.1, 
for which our sample is spatially complete on scales of 
23.7 < R prop < 29.7 h^ 1 kpc. In Figure [4] we compare 
the small-scale clustering of quasars determined from our 
complete PR4 binary quasar sample to the results from 
H06. In the left-hand panel, we consider our calcula- 
tion for all quasars in Table [H and project the result for 
both proper and comoving scales by placing all binaries 
at the mean proper or comoving distance of our sample. 
In the right-hand panel, we consider binaries in Table [4] 
with 23.7 < R pr0 p < 29.7 h" 1 kpc and 1.0 < z < 2.1. 
We additionally include SPSSJ1635+2911 in this second 
sample, as we are spatially complete for component sep- 
arations of 23.7 < Rprop < 29.9 for 1.03 < z < 2.1. 

Figure [5] demonstrates that our results, using a com- 
plete, statistical sample of UVX binary quasars are 
broadly consistent with H06. To compare with H06 we 
fit power laws, displayed in Figure [H out to proper (co- 
moving) scales of < 100 h~ x kpc (< 200 h^ 1 kpc), and 
integrate them over our scales of interest. At the mean 
redshift (z = 1.40) of our sample, 3.9" < 9 < 5.2", the 
range for which our sample is complete, is equivalent to 
proper (comoving) scales of 23.5 h~ x kpc < R prop < 
31.4 h- 1 kpc (56.5 h- 1 kpc < R < 75.3 h^ 1 kpc). Over 
this angular range, we find W p — 30.8±}f p for our data. 
The proper (comoving) power-law fit to the data from 
H06, implies W p = 55.1 (60.4) over the same scales, a 
1.5<t (1.9cr) difference. The difference is slightly more 
pronounced if we determine W p for the H06 data at the 
mean scale of the 8 binaries in our sample, instead of pro- 
jecting 3.9" < 6 < 5.2" back to z = 1.4. For the 5 bina- 
ries in our spatially complete clustering subsample, which 
covers scales of 23.7 hr 1 kpc < R pr0 p < 29.9 hr 1 kpc, 
we find W p = 24.0±\fl, and the H06 data implies 
Wp = 57.2, 2.4 times (and 2.0c) higher than our result. 

5. DISCUSSION 

We find that the projected correlation function of 
quasars at proper scales of ~ 25 h~ x kpc has an am- 
plitude a factor of 2.4 times lower than that deter- 
mined by H06. H06 argue that clustering on these scales 
is ~ 10 times higher than e xpected from projec ting 
the quasar autocorrelation of iPorciani et al.l (|2004h to 
smaller scales. Our data thus imply excess quasar clus- 
tering at ~ 25 h^ 1 kpc of a factor of ~ 4, consistent with 
the quoted upper limit in M07b. Given that our sample 
is targeted differently than any previous samples of bi- 
nary quasars (i.e. UVX but otherwise regardless of the 
color similarity of the candidate components) , and given 
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our simple selection function, our work might be viewed 
as independently corroborating the evidence for excess 
quasar clustering on small scales first detected in H06. 
We note that our binary quasar clustering subsample is 
largely independent of the sample used by H06, as of the 
8 binaries listed in Tabled] only 2 appear in Table [21 

we find a significantly smaller 
clustering amplitude than found by H06. As the bina- 
ries in Table 2] are all at z < 2 but H06 include sev- 
eral Rprop ~ 25 hr 1 kpc binaries with z > 2, an in- 
teresting possibility is that W p is a function of scale 
and redshift. Certainly, on large scales (> 1 hr 1 Mpc), 
quasa rs cluster twice as strongly at z > 2 t han at z < 1 
(e.g.. lCroom~et"aH [20051: IShen et al J 12001 . The mean 
redshift of the H06 sample is I = 1.87 (J. Hennawi 
2007, private communication ) compared t o z = 1.42 
for our sample in Table [U ICroom et al.l (|2005l ) esti- 
mate quasar bias on large (i.e., > 1 h^ 1 Mpc) scales 
to follow b(z) = 0.53 + 0.289(1 + zf . Given that quasar 
clustering scales as b 2 , the implied relative amplitude be- 
tween quasar clustering in H06 and our sample is ~ 1.72. 
Scaling our measurement of W p — 30.8±}f by this fac- 
tor implies W p — 53.0±f| - g, easily consistent with the 
value of W p — 57.2 implied by our fit to H06 in Fig- 
ure 2] Scaling the W p — 24.0±}oJ estimate for our 
complete subsample (z — 1.60) in the same way implies 
W p — 33.1±i4;|, lower than W p = 57.2 but consistent 
to ~ la. Thus, multiplying our quasar clustering am- 
plitudes on small scales by the implied bias evolution 
from quasar clustering on large scales somewhat recon- 
ciles our clustering amplitudes with the higher redshift 
results from H06. Consistent evolution of quasar clus- 
tering on large and small scales is a natural feature of 
merger-driven models of quasar origin (e.g., Figure 17 of 
iHopkins et aTll2007bl ). 

We can also look for direct redshift dependence within 
our sample. Splitting the sample in Table[4]at z — 1.55, 
wc find Wp = 27.9±f|;| for z < 1.55 {z ~ 1.07) and 
Wp = 34.9±?y for z > 1.55 (z ~ 1.76). Although these 
results seem consistent, the mean separation of the bina- 
ries at z < 1.55 (z > 1.55) is Rp rop = 23.8 (28.4). The 
power law fits displayed in Figure 3] imply that we should 
therefore expect quasars to be ~ 30% more clustered at 
z < 1.55 than at z > 1.55. Incorporating this ~ 30% 
correction, we measure the ratio of the amplitudes of W p 
at z > 1.55 and z < 1.55 to be 1.62 ± 1.54. Although 
not a very significant detection, a factor of 1.62 ± 1.54 
is again consistent with quasar clustering evolution on 
large scales, as a b(z) — 0.53 + 0.289(1 + z) 2 bias relation 
implies a ratio of ~ 2.4 between the amplitudes of W p at 
z = 1.76 and z = 1.Q7. 

Our analyses could be easily extended, as the 98 
(3" < AO < 6") DR4 KDE binary quasar candidates 
should contain ~ 40 binary quasars. Based on our W p 
estimates, such a sample is the minimum necessary to 
detect any redshift dependence to quasar clustering at 
Rprop ~ 25 h~ x kpc, providing only a ~ 1.5er detec- 
tion. Completing this sample on a 4-meter class tele- 
scope would likely require dark time, as 14 of the DR4 
binar y candidates have a component at g > 21. Alterna- 
tively, IShen et all (|2007f) suggest a survey of quasar pairs 
at z > 3, where stronger quasar clustering may lead to a 



more significant detection of evolution. The DR6 KDE 
quasar catalog (Richards et al, in prep), which will be 
selected to i < 21 (g < 21.25), should contain ~ 500 
binary quasar candidates across a large redshift range 
(0.4 < z < 5) with which to pursue this goal in a single 
homogeneous sample. 

6. CONCLUSIONS 

We have presented a spectroscopic survey of binary 
quasar candidates with separations in the range 3" < 
AO < 6". Our candidates (see Table 1 of M07b) are a 
subsample of quasars photometri cally classified in SPS S 
DR4 using the KDE technique of lRichards et all (|2004h . 
We define a binary as a quasar pair with a line-of-sight 
velocity difference of |Aun| < 2000 km s _1 (see H06). 
We present 27 new sets of observations and identify both 
members of 22 candidate binary quasars. Of the 22 new 
pairs, 10 turn out to be binary quasars (of which ~ 2 
might actually be lenses) . This roughly doubles the num- 
ber of known binary quasars with 3" < A8 < 6" at 
z < 2 and g < 21. A further 9 of our observed candi- 
dates are projected quasar pairs, and 3 contain a NELG 
or star. This confirms that the KD E technique is ~ 95% 
efficient at selecting quas ars (e.g., iRichards et al.l [20041 : 
iMvers et al.l 120061 l2007af) . Combined with observations 
from the literature (mainly from H06), 46% of the DR4 
KDE binary quasar candidates have now been observed, 
of which ~ 47% are binaries or lenses, ~ 40% are pro- 
jected quasar pairs, and the remainder contain a non- 
quasar. 

As our candidate binaries, beyond a UVX cut, are se- 
lected regardless of the relative colors of the quasars in 
the pair, we can try to assess the color similarity criteria 
that optimally select binary quasars. For quasars pre- 
selected with an already efficient approach such as the 
KDE technique, we find that a % 2 color similarity statis- 
tic of xloior < 70 wm return 97% of binaries, lenses and 
projected pairs (multiplied by the 95% completeness of 
the KDE technique itself) at 93% efficiency. Most of this 
efficiency in selecting quasar pairs comes from the KDE 
technique itself, as imposing no color similarity criterion 
is 87% efficient. To select binary quasars, while reject- 
ing projected quasar pairs, we recommend a cut in the 
overlap of the photometric redshifts of the two candidate 
quasars in a pair. An overlap of > 0.03 in the primary 
solution for the photometric redshift probability density 
functions of the pair can be constructed to be ~ 95% 
complete and ~ 70% efficient for binary quasars. Sim- 
ilarly, of course, the reverse probability cut of < 0.03, 
perhaps coupled with a Xcolor < 70 cut to remove stars, 
can be used to reject binary quasars in favor of projected 
quasar pairs. 

We measure the clustering of a complete sample of DR4 
binaries on proper scales of 23.7 < R pr0 p < 29.7 kpc. 
We find that, at ~ 25 h~ l kpc, quasars cluster with an 
amplitude 2.4 times, or 2. Oct, lower than determined by 
H06. As the mean redshift of the H06 sample is z — 1.87 
compared to z — 1.40 for our sample, this can be inter- 
preted as evidence of evolution in quasar clustering on 
scales of ~ 25 hr 1 kpc. The implied evolution is broadly 
consistent with merger-driven models, where the quasar 
population is expected to ev olve with consistent large- 
to-small scale clustering (e.g., [Hopkins et~aT1l2007bl) . We 
find no significant evidence for quasar clustering evolu- 
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tion at ~ 25 h" 1 kpc, from z ~ 1.07 to z ~ 1.76, in our 
sample alone. Assuming evolution in the binary quasar 
population at the level suggested by our current sample, 
we argue that the 98 quasars in the DR4 KDE candi- 
date binary sample will detect any clustering evolution 
at proper scales of ~ 25 hr 1 kpc at ~ 1.5c significance. 
A sample of ~ 200 candidates (~ 80 binary quasars) will 
be necessary to definitively detect clustering evolution at 
~ 25 h' 1 kpc for z < 2.5. 
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Fig. 1. — Two of the five quasar pairs in our DR4 KDE sample that require higher resolution spectroscopy to determine whether they are 
binary quasars or a lensed quasar. These spectra were taken with the R-C Spectrograph on the Mayall 4-m at KPNO at a resolution ~ 5A, 
and have been smoothed with a 5 pixel boxcar. SDSSJ1158+1235A,B and SDSSJ1606+2900A,B are our two most likely lens candidates, 
based on their similar colors and component separations of AS < 4". In both panels, component A is the upper spectrum and component 
B has been offset, which can cause component B to have a flux density below zero. In each case, both components are, in reality, at nearly 
identical flux levels. 
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Fig. 2.— Spectra of SDSSJ1235+6836A,B taken with the R-C Spectrograph on the Mayall 4-m at KPNO at a resolution ~ 5A, and 
smoothed with a 5 pixel boxcar. In the lower panel, the vertical lines mark common quasar emission lines. For component A, in the upper 
panel, emission lines are marked at the systemic redshift of com ponent B offset by ±2000 km s — 1 , a typical window for a binary quasar 
(e.g., H06). The systemic redshift of component A in Table ITdl is derived assuming the emission line near 3900 A is CIV, and is close to 
the red side of these windows (z = 1.529). At red wavelengths, component A resembles an M star. 
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Fig. 3. — Upper-left: The cumulative probability distribution of the x\ i or c °l° r similarity statistic (Equation[TJ for the 45 spectroscopically 
observed DR4 KDE candidate binary quasars, the 53 that are not yet observed, and for all 98 candidate pairs. Upper-right: Similar to 
upper-left, but for the 21/45 observed candidates that are binary quasars (or lenses), for the 24/45 that are not binaries and for the 18/45 
quasar pairs at disjoint redshifts. Lower-left: The selection completeness and efficiency effects of imposing a X^oior urQ it on * ne 45 observed 
DR4 KDE candidate binaries. Efficiency is (number of binaries < Xcoior)/(*°'' a ^ candidates < xtnlor)' completeness is (number of binaries 
or lenses < X 2 i„-.)/21. Lower-right: Similar to lower-left, but using the fractional overlap of the photometric redshift solutions for the 
candidate components (z p hot) as the determinant of binary selection. Note that more overlap in photometric redshift implies greater color 
similarity. 
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Fig. 4. — The projected correlation function, W p , of quasars on small scales. Left: The results of H06 in proper (lower axis) and comoving 
(upper axis) coordinates are compared to our mean DR4 KDE spectroscopic results averaged over 3.9" < A9 < 5.2" and projected to 
a mean transverse separation at z = 1.4 (the mean redshift of both the DR4 KDE quasar sample and the subset of 16 quasars listed in 
Table U), in proper (open triangle) and comoving (solid triangle) coordinates. The dashed (solid) line shows a power law fit to the H06 
data on proper (comoving) scales < 100 h~ 1 kpc (< 200 h —1 kpc). Right: A similar comparison but our data are now averaged over scales 
Rvrop < 29.7 h- 1 kpc and redshifts of 1.0 < z < 2.1. The DR4 KDE binary quasar sample is spectroscopically complete for 
< 20.85, and 3.9" < Ad < 5.2", and for 23.7 < Rprop < 29.7 h" 1 kpc over 1.0 < z < 2.1. 
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TABLE 1a 

DR4 KDE CANDIDATE BINARIES FOR WHICH WE HAVE OBSERVED ONE MEMBER 



A6» 


y 2 


Name 


a (J2000) 


S (J2000) 


g 


z phot range 


SDSS z 


Our z 


5.50 


36.2 


SDSSJ093014.81+420038.7 
SDSSJ093015.01+420033.6 


09 30 14.816 
09 30 15.016 


+42 00 38.71 
+42 00 33.68 


19.71 
20.03 


1.15,1.425,1.50,0.93 
1.85,2.025,2.15,0.62 




0.544 7 


5.55 


24.2 


SDSSJ093521. 02+641219. 8 
SDSSJ093521. 80+641221. 9 


09 35 21.020 
09 35 21.807 


+64 12 19.89 
+64 12 22.00 


20.99 
20.96 


1.45,1.775,2.10,0.92 
0.25,0.375,0.45,0.66 




1.566 


5.33 


6.4 


SDSSJ095840.74+332216.3 
SDSSJ095840.94+332211.5 


09 58 40.746 
09 58 40.945 


+33 22 16.31 
+33 22 11.59 


19.18 
20.64 


1.35,1.475,2.10,0.84 
1.45,1.725,2.10,0.91 


1.891 


1.888 


3.46 


161 


SDSSJ103939.31+100253.0 
SDSSJ103939.53+100254.3 


10 39 39.317 
10 39 39.532 


+ 10 02 53.01 
+ 10 02 54.40 


18.42 
19.60 


0.10,0.175,0.25,0.56 
0.65,1.125,1.55,0.92 


0.161 


0.161 


4.22 


8.2 


SDSSJ162847.75+413045.4 
SDSS J162848.06+413043. 1 


16 28 47.752 
16 28 48.069 


+41 30 45.45 
+41 30 43.19 


19.81 
20.40 


1.35,1.525,1.70,0.85 
1.95,2.075,2.20,0.63 




0.831 ? 




Note 


. — The angular pair separations arc denoted A8 ( ), 


and X^olor * s eacn 


pair's 


color similarity statistic 


(Equation [T}. The 



photometric rcdshifts z p hot are expressed as "lowest extent, peak, highest extent, probability of true redshift lying in this range". The 
"SDSS z" column shows matches to any spectroscopic object in the SDSS DR6 Catalog Archive Server (mainly, e.g., Schneider ct al. 
2007). The "Our z" column lists our new spectroscopic confirmations from KPNO data. In the "SDSS z" and "Our z" columns the 
object is a quasar at the provided redshift, unless otherwise noted. Rcdshifts labeled with a ? arc based on a single emission line, 
which is reasonably assumed to be Mgll. g is not corrected for Galactic extinction. 



TABLE 1b 

Confirmed binary quasars in the DR4 KDE candidate sample 



A8 


Y 2 , 


Name 


a (J2000) 


8 (J2000) 


9 


Zphot range 


SDSS z 


Our z 


*3.56 


4.0 


SDSSJ1158+1235A 
SDSSJ1158+1235B 


11 58 22.776 
11 58 22.989 


+12 35 18.59 
+12 35 20.31 


19.90 
20.12 


0.45,0.525,0.65,0.55 
0.40,0.475,0.65,0.63 




0.596 7 
0.596 ? 


*4.74 


32.0 


SDSSJ1320+3056A 
SDSSJ1320+3056B 


13 20 22.545 
13 20 22.643 


+30 56 22.87 
+30 56 18.29 


18.60 
19.92 


1.35,1.475,1.65,0.90 
1.45,1.575,1.80,0.51 


1.597 


1.595 
1.596 


*4.50 


2.8 


SDSSJ1418+2441A 
SDSSJ1418+2441B 


14 18 55.418 
14 18 55.536 


+24 41 08.92 
+24 41 04.71 


19.27 
20.22 


0.45,0.525,0.70,0.96 
0.40,0.625,0.70,0.86 


0.573 


0.572 
0.573 


4.27 


3.2 


SDSSJ1426+0719A 
SDSSJ1426+0719B 


14 26 04.266 
14 26 04.326 


+07 19 25.86 
+07 19 30.04 


20.82 
20.12 


0.95,1.175,1.45,0.99 
1.00,1.225,1.45,0.97 




1.312 
1.309 


5.41 


25.6 


SD9SJ1430+0714A 
SDSSJ1430+0714B 


14 30 02.664 
14 30 02.886 


+07 14 15.62 
+07 14 11.33 


20.27 
19.50 


1.00,1.225,1.40,0.97 
1.05,1.375,1.45,0.97 


1.258 


1.246 
1.261 


5.14 


65.6 


SDSSJ1458+5448A 
SDSSJ1458+5448B 


14 58 26.165 
14 58 26.728 


+54 48 14.85 
+54 48 13.19 


20.79 
20.53 


1.50,1.775,1.95,0.75 
1.65,1.925,1.95,0.47 




1.913 
1.912 


"3.45 


10.8 


SDSSJ1606+2900A 
SDSSJ1606+2900B 


16 06 02.812 
16 06 03.021 


+29 00 48.79 
+29 00 50.88 


18.50 
18.42 


0.50,0.725,1.00,0.65 
0.70,0.875,1.00,0.92 


0.770 


0.769? 
0.769? 


*4.92 


25.6 


SDSSJ1635+2911A 
SDSSfJ1635+2911B 


16 35 10.148 
16 35 10.306 


+29 11 20.65 
+29 11 16.19 


18.83 
20.43 


1.45,1.575,1.80,0.84 
1.40,1.525,1.85,0.79 


1.586 


1.582 
1.590 





NOTE. — We define a binary quasar by a line-of-sight velocity difference of |Av||[ < 2000 km s 1 in the rest-frame of 
either component. The pair containing SDSS J143002. 66+071415. 6 has velocity difference |At>||| — 2000 ± 400 km s~ 1 , just 
inside our definition of a binary. Components arc denoted A and B so that the position angle from A to B lies between 0° 
and 180°. A preceding * denotes that our spectroscopy alone is insufficient to r ule out a lens interpretati on for this pair (see 
32.2.11 . SDSSJ1320+3056A first appeared with a confirmed redshift (z=1.587) in lVeron-Cettv et all d2004fl . g is not corrected 
for Galactic extinction. See Table [Ta\ for additional notes describing shared notation. 
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TABLE lc 

Confirmed projected pairs in the DR4 KDE candidate sample 



A8 Xcolor Name a (J 2000 ) S (J2000) g z phot range SDSS z Our z 



4.28 1220 SDSSJ083258.34+323003.3 08 32 58.348 +32 30 03.35 19.96 2.75,2.775,2.80,0.88 0.397 

SDSSJ083258.56+323000.0 08 32 58.567 +32 30 00.08 19.59 0.40,0.425,0.50,0.70 star 

3.42 2.0 SDSSJ084257.37+473342.5 08 42 57.378 +47 33 42.56 19.00 0.50,0.625,0.70,0.51 1.552 1.554 

SDSSJ084257.63+473344.7 08 42 57.638 +47 33 44.74 20.45 0.85,1.775,2.10,0.84 1.681 

4.28 4.9 SDSSJ085914.77+424123.6 08 59 14.771 +42 41 23.67 21.02 0.95,1.375,1.65,0.66 1.396 

SDSSJ085915.15+424123.5 08 59 15.159 +42 41 23.58 19.22 0.80,0.975,1.40,0.93 0.898 0.902 ? 

X57 206 SDSSJ094309.36+103401.3 09 43 09.363 +10 34 01.31 20(37 0.95,1.525,1.70,0.80 1.431 

SDSSJ094309.66+103400.6 09 43 09.670 +10 34 00.65 19.24 1.05,1.325,1.40,0.99 1.238 1.240 

"5T91 3T8 SDSSJ111053.63+605347.9 11 10 53.633 +60 53 47.97 TSM 0.65,0.775,0.90,0.96 0.793 

SDSSJ111054.10+605343.1 11 10 54.105 +60 53 43.16 20.86 1.55,1.775,2.20,0.68 0.552 7 

5.19 18.4 SDSSJ112556.32+143148.0 11 25 56.321 +14 31 48.10 20.70 0.90,1.325,1.60,0.69 NELG 

SDSSJ112556.54+143152.1 11 25 56.549 +14 31 52.10 20.38 1.50,1.725,2.15,0.97 1.924 

5.78 29.9 SDSSJ113637.52+563500.4 11 36 37.526 +56 35 00.48 19.94 1.85,2.075,2.15,0.63 1.282 

SDSSJ113638.09+563503.9 11 36 38.090 +56 35 03.90 19.02 0.60,0.675,0.90,0.56 2.672 

4.92 39.3 SDSSJ114503.06+660211.3 11 45 03.063 +66 02 11.35 20.09 1.65,1.825,2.00,0.89 1.732 

SDSSJ114503.74+660208.6 11 45 03.741 +66 02 08.67 20.24 0.50,0.675,0.85,0.52 2.304 

"3\26 1TT5 SDSSJ123122.27+493433.8 12 31 22.276 +49 34 33.86 20.63 0.50,0.725,0.90,0.53 0.780' 

SDSSJ123122.37+493430.7 12 31 22.378 +49 34 30.75 19.94 1.55,1.775,2.15,0.95 1.811 ? 

TOO 107 SDSSJ150656.86+505610.5 15 06 56.866 +50 56 10.56 19.22 0.65,0.775,1.00,0.75 0.775' 
SDSSJ150657.18+505607.9 15 06 57.183 +50 56 07.92 19.75 2.00,2.225,2.40,0.45 2.204 



Note. — SDSSJ112556. 32+ 143148.0, the NELG, has redshift z=0.246. Subsequent to our observations, SDSSJ094309. 36+103401. 3 
appeared in llnada et a.1.1 J2007t> . with z = 1.433. g is not corrected for Galactic extinction. The redshift for SDSS J123122. 37+493430. 7 
is based on a single CIV emission line, with weak confirming CIII], See Table [Ta] for additional notes describing shared notation. 



TABLE Id 

Ambiguous pairs in the DR4 KDE binary quasar candidate sample 



A8 X 2 co i or Name a (J2000) <5 (J2000) g z phot range SDSS z Our z 



4.80 


42.5 


SDSSJ093424.11+421135.0 
SDSSJ093424.32+421130.8 


09 34 24.112 
09 34 24.324 


+42 11 35.05 
+42 11 30.87 


21.01 
20.30 


1.45,1.675,2.20,0.90 
1.00,1.125,1.40,0.99 


f/lcss 
1.339 


3.95 


10.7 


SDSSJ120727.09+140817.1 
SDSSJ120727.25+140820.3 


12 07 27.100 
12 07 27.259 


+ 14 08 17.18 
+ 14 08 20.38 


20.39 
20.27 


1.60,1.775,2.00,0.89 
1.55,1.775,1.95,0.86 


1.801' 

1.8" 


3.51 


207 


SDSSJ1235+6836A 
SDSSJ1235+6836B 


12 35 54.783 
12 35 55.270 


+68 36 24.78 
+68 36 27.07 


19.04 
19.70 


2.75,2.775,2.80,0.55 
0.50,0.625,1.10,0.51 


1.529" 
1.514 


4.35 


14.4 


SDSSJ1507+2903A 
SDSSJ1507+2903B 


15 07 46.909 
15 07 47.234 


+29 03 34.15 
+29 03 33.28 


20.44 
19.97 


0.80,0.975,1.25,0.77 
0.70,0.775,0.95,0.66 


0.875- 
0.862 ? 





Note. — Redshifts marked ?? are derived from a single emission line. This differs from the ? notation as the redshift is based 
on similar emission in the other component (rather than simply assuming that the emission is Mgll). The ambiguities, and why 
we conclude that SDSSJ1235+6836 and SDSSJ1507+2903 are binaries but the other pairs are not, arc discussed in i|2.2.2| g is not 
corrected for Galactic extinction. See Table \Ta\ for additional notes describing shared notation. 
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TABLE 2 

Previously identified DR4 KDE binary quasar candidates (3" < Ad < 6") 



A0 


* 2 , 


Name 


a (J 2000) 


o (J 2000) 


9 


z P hot range 


z 


Projected pairs 
4.90 36.7 


SDSSJ024907.77+003917.1 
SDSSJ024907.86+003912.4 


02 49 07.778 
02 49 07.866 


+00 39 17.12 
+00 39 12.40 


19.36 
20.63 


2.00,2.175,2.25,0.48 
0.45,0.675,0.85,0.84 


2.164 
star 


4.09 


18.0 


SDSSJ083649. 45+484150.0 
SDSSJ083649.55+484154.0 


08 36 49.456 
08 36 49.554 


+48 41 50.08 
+48 41 54.06 


19.31 
18.50 


0.45,0.675,0.80,0.67 
1.50,1.675,1.95,0.94 


0.657 
1.712 


5.42 


4.5 


SDSSJ090235. 35+563751. 8 
SDSSJ090235.73+563756.2 


09 02 35.356 
09 02 35.731 


+56 37 51.84 
+56 37 56.29 


20.95 
20.56 


1.15,1.275,1.45,0.98 
1.05,1.225,1.45,0.98 


1.39 
1.34 


3.14 


37.1 


SDSSJ095454.73+373419.7 
SDSSJ095454.99+373419.9 


09 54 54.735 
09 54 54.999 


+37 34 19.79 
+37 34 19.99 


19.57 
18.91 


0.95,1.475,1.65,0.90 
1.45,1.575,1.95,0.94 


1.544 
1.892 


4.76 


10.2 


SDSSJ114718.44+123439.8 
SDSSJ114718.66+123436.3 


11 47 18.448 
11 47 18.668 


+12 34 39.84 
+12 34 36.33 


20.91 
19.80 


1.45,1.625,2.00,0.66 
2.15,2.225,2.60,0.54 


1.583 
2.232 


3.06 


5.6 


SDSSJ120450.54+442835.8 
SDSSJ120450.78+442834.2 


12 04 50.543 
12 04 50.784 


+44 28 35.89 
+44 28 34.25 


19.04 
19.48 


0.95,1.125,1.45,0.98 
1.35,1.725,1.95,0.78 


1.144 
1.814 


5.04 


23.3 


SDSSJ124948.12+060709.0 
SDSSJ124948. 17+060714.0 


12 49 48.127 
12 49 48.179 


+06 07 09.04 
+06 07 14.02 


20.41 
20.38 


2.20,2.325,2.65,0.89 
1.85,2.075,2.20,0.58 


2.001 
2.376 


3.01 


20.0 


SDSSJ125530.44+630900.5 
SDSSJ125530.82+630902.0 


12 55 30.445 
12 55 30.823 


+63 09 00.51 
+63 09 02.09 


20.30 
20.60 


1.50,1.675,1.90,0.88 
1.10,1.375,1.50,0.98 


1.753 
1.393 


4.94 


31.2 


SDSSJ142359.48+545250.8 
SDSSJ142400.00+545248.7 


14 23 59.484 
14 24 00.006 


+54 52 50.83 
+54 52 48.79 


18.63 
19.93 


1.00,1.175,1.45,0.973 
1.45,1.575,1.90,0.772 


1.409 1 
0.610 1 


4.35 


25.2 


SDSSJ162902.59+372430.8 
SDSS J162902.63+372435. 1 


16 29 02.594 
16 29 02.634 


+37 24 30.85 
+37 24 35.17 


19.17 
19.35 


0.80,0.975,1.10,0.93 
0.70,0.925,1.10,0.98 


0.923 
0.906 


5.83 


176 


SDSSJ171334.41+553050.3 
SDSSJ171335.03+553047.9 


17 13 34.414 
17 13 35.037 


+55 30 50.36 
+55 30 47.91 


18.88 
19.11 


1.00,1.375,1.45,0.968 
2.00,2.175,2.20,0.686 


1.276 1 
star 1 


5.06 


172 


SDSSJ211157.24+091559.3 
SDSSJ211157.26+091554.2 


21 11 57.248 
21 11 57.269 


+09 15 59.33 
+09 15 54.28 


20.73 
19.83 


0.95,1.275,1.40,0.995 
1.00,1.325,1.35,0.983 


star s 


Binary quasars 
3.94 18.8 


SDSSJ0959+5449A 
SDSSJ0959+5449B 


09 59 07.060 
09 59 07.471 


+54 49 08.09 
+54 49 06.38 


20.60 
20.07 


1.90,2.025,2.15,0.59 
1.40,1.575,2.10,0.90 


1.956 
1.954 


3.55 


6.5 


SDSSJ1259+1241A 
SDSSJ1259+1241B 


12 59 55.464 
12 59 55.617 


+12 41 51.06 
+12 41 53.81 


19.99 
20.09 


1.95,2.175,2.30,0.43 
1.90,2.175,2.25,0.52 


2.180 
2.189 


3.81 


4.3 


SDSSJ1303+5100A 
SDSSJ1303+5100B 


13 03 26.144 
13 03 26.177 


+51 00 51.00 
+51 00 47.21 


20.54 
20.37 


1.50,2.075,2.20,0.82 
1.60,1.775,2.00,0.93 


1.686 
1.684 


3.12 


0.5 


SD9SJ1337+6012A 
SDSSJ1337+6012B 


13 37 13.085 
13 37 13.133 


+60 12 09.70 
+60 12 06.60 


20.04 
18.59 


1.30,1.775,2.05,0.66 
1.50,1.625,1.95,0.94 


1.721 
1.727 


5.13 


4.1 


SDSSJ1432-0106A 
SDSSJ1432-0106B 


14 32 28.949 
14 32 29.247 


-01 06 13.55 
-01 06 16.06 


21.10 
17.83 


1.55,2.125,2.25,0.69 
1.90,2.025,2.15,0.96 


2.082 
2.082 


4.11 


7.2 


SDSSJ1530+5304A 
SDSSJ1530+5304B 


15 30 38.564 
15 30 38.824 


+53 04 04.03 
+53 04 00.65 


20.56 
20.70 


1.45,1.575,1.95,0.65 
1.40,1.725,2.15,0.93 


1.531 
1.533 


3.90 


2.8 


SDSSJ1637+2636A 
SDSSJ1637+2636B 


16 37 00.881 
16 37 00.932 


+26 36 13.71 
+26 36 09.87 


20.61 
19.36 


0.45,0.575,0.85,0.46 
1.40,1.525,1.80,0.64 


1.961" 
1.961° 


3.72 


3.9 


SDSSJ1723+5904A 
SDSSJ1723+5904B 


17 23 17.307 
17 23 17.421 


+59 04 42.79 
+59 04 46.41 


20.31 
18.88 


1.45,1.725,2.25,0.63 
1.55,1.725,1.90,0.94 


1.597 
1.604 


5.81 


35.0 


SDSSJ2214+1326A 
SDSSJ2214+1326B 


22 14 26.792 
22 14 27.032 


+13 26 52.38 
+13 26 57.01 


20.64 
20.34 


1.55,2.025,2.20,0.85 
1.65,1.825,2.05,0.96 


1.995 
2.002 


Confirmed lenses 

3.76 1.9 SDSSJ1004+4112A 
SDSSJ1004+4112B 


10 04 34.800 
10 04 34.917 


+41 12 39.29 
+41 12 42.81 


18.64 
19.04 


1.55,1.725,2.00,0.94 
1.55,1.725,2.15,0.79 


1.734 1 
1.734 1 


3.04 


46.4 


SDSSJ1206+4332A 
SDSSJ1206+4332B 


12 06 29.648 
12 06 29.652 


+43 32 17.57 
+43 32 20.61 


18.78 
19.38 


1.65,1.825,2.05,0.96 
1.95,2.175,2.35,0.53 


1.789° 
1.789° 





NOTE. — Components of a binary are denoted A and B so that the position angle from A to B lies between 0° and 180°. 
This convention differs from H06, from which we take identifications and rcdshifts, e xcept for objects labeled S (t aken from 
the SDSS), D (discove red by [SraJne^fc^^^dman| 1978|, confirmed as a poss ible lens by Diorgovski & Spinrad 198 4|, and likel y 
a binary instead, e.g . jKocha^ e^^^^![[^9 ^r[Peng ct al. 1999; Rusinl 120021') . i (part of the quad lens from llnada et al.ll2003l) , 
o l lOguri et all 120051) , and I <Inada et alj|2007l) . Both quasars SDSSJ162902 59+37243 .8 and SDSS J162902. 63+372435. 1 
first appear in IMason et al.ll2000D . SDSSJ1004+4112A was discovered by ICao et al] fll999D , and SDSS J1432-0106B by 
IHewett et al.l J1991I) . We note that we mistakenly listed SDSS J095454. 73+373419. 7 as lying at z=1.554 in M07b. g is not 
corrected for Galactic extinction. See Table [Ta] for additional notes describing shared notation. 
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TABLE 3 

Breakdown of (3" < A6» < 6") DR4 KDE quasar pairs 



Category Number of Confirmed Pairs 

Total binary quasar candidates 98 

Total now identified 45 

Likely binary quasars 19 

Quasar pairs separated in redshift 18 

Pairs containing > 1 non-quasars 6 

Pairs that are confirmed lenses 2 



Note. — It is possible that a few more objects listed as "Likely 
binary quasars" may turn out to be a lensed quasar when scrutinized 
at higher resolution. 



TABLE 4 

Complete, statistical, clustering subsample 



Name A0(") xLor R P r °P R Z A Z B |Al)||| Table 

SDSSJ0959+5449 
SDSSJ1320+3056 
SDSSJ1418+2441 
SDSSJ1426+0719 
SDSSJ1458+5448 
SDSSJ1507+2903 
SDSSJ1530+5304 
SDSSJ1635+2911 




Note. — The DR4 KDE binary quasar candidate sample is now spectroscopically complete 
for component separations 3.9" < A6 < 5.2" for g < 20.85 in regions with Galactic absorp- 
tion A g < 0.17. 20 pairs meet these criteria, and 8 of them are (the listed) binary quasars. 
A * denotes a possible lens (see note in Table [Tbl . R prov (R) is the transverse proper (co- 
moving) separation (h" 1 kpc). | Ally | is the line-of-sight velocity difference (km s _1 ). The 
final column "Tabic" denotes where we first listed these binaries. The 5 listed binaries with 
transverse separations of 23.7 < R prop < 29.9 represent a spatially complete subsample for 
redshifts of 1.03 < z < 2.10. 



